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ABSTRACT 


We  report  one-  (OP A)  and  two-photon  absorption  (TP A)  excitation  energies  and  cross  sections 
for  a  series  of  7-aminocoumarins  using  time-dependent  density  functional  theory  with  various  basis  sets 
and  functionals,  including  exchange-correlation  functionals  using  the  Coulomb-attenuating  method 
(CAMB3LYP  and  mCAMB3LYP),  to  evaluate  their  performance  in  the  gas -phase  and  in  solvents. 
Except  for  the  CAMB3LYP  results,  the  computed  OPA  excitation  energies  and  transition  dipole 
moments  are  in  good  agreement  with  experiment.  The  range  of  errors  obtained  from  various 
functionals  is  discussed  in  detail.  The  relationship  of  donor  and  acceptor  groups  with  the  OPA  and  TPA 
resonances  and  intensities  is  also  discussed. 
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1.  Introduction 


Coumarin  derivatives  that  are  formed  by  various  substitutions  of  amino  groups  at  the  7-position  are 
widely  used  as  probes  in  ultrafast  solvation  studies,  as  laser  dyes,12  and  as  nonlinear  optical 
chromophores3'5  in  optoelectronic  devices.  Thus,  pertinent  optical  properties  for  these  applications  have 
been  extensively  reported  in  the  literature.  Electronic  one-photon  absorption  (OPA)  spectra  of 
coumarins,  in  particular,  have  been  characterized  in  high  temperature  vapors,6'  in  a  supersonic  jet 
expansion,6"11  and  in  a  variety  of  polar  and  nonpolar  solvents.9'1'1  Experimental  studies  have  revealed  the 
polar  nature  of  the  first  excited  state,  which  has  been  characterized  as  an  intramolcular  charge-transfer 
(ICT)  state  with  a  larger  dipole  moment  than  in  the  ground  state.  The  large  difference  in  dipole  moment 
between  the  ground  and  excited  states  (A/ige)  provides  a  viable  mechanism  for  the  intense  two-photon 

absorption  (TPA)16'19  and  other  nonlinear  processes3  4  that  have  been  observed  and  quantified  for  a 
number  of  coumarin  derivatives. 

Theoretically,  OPA  properties  of  coumarins  in  the  gas-phase  have  been  reported  in 
semiemperical  studies  using  AMI  ,  PM3  ’  and  MNDO  ‘  methods.  Recently,  Cave  et  al.  reported  a 
detailed  theoretical  study  of  coumarin  120  (C120)  and  C151  using  a  wide  range  of  electronic  structure 
methods  to  evaluate  various  ground  and  excited  state  properties.  They  found  that  using  time -dependent 
density  functional  theory  (TDDFT)25'28  with  hybrid  functionals  gives  excellent  agreement  with  So-Si 
experimental  OPA  excitation  energies.  In  another  survey29  involving  Cl 02,  Cl 52,  Cl 53,  and  C343, 
excellent  agreement  with  experiments  was  also  reported  for  TDDFT  excitation  energies.  Optical 
properties  of  Cl  20, 24  C151,24’30'32  1  52A  (C35,  C481)30*31  and  C15330’31  have  been  studied  using  TDDFT 
with  continuum,  and  explicit"  solvent  models.  More  recently,  Ingrosso  et  al.  reported  "  a  solvation 
study  of  Cl 53  in  supercritical  fluoroform  using  TDDFT  with  Becke's  three-parameter  (B3LYP)’l  ,h 
hybrid  functional.  The  successful  applications  of  TDDFT  with  hybrid  functionals  for  coumarins  are 
consistent  with  our  previous  studies'’7  on  the  prediction  of  excitation  energies  involving  charge-transfer 
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excited  states.  However,  we  found  the  B3LYP  functional34"36  and  the  hybrid  functional  of  Perdew, 
Burke,  and  Erzerhof  (PBEO  underestimate  the  experimental  excitation  energies  of  a  number  of  ICT 
chromophores.  Both  functionals  were  found  to  overestimate  the  oscillator  strengths  for  the  first  ICT 
excited-states.  The  failures  to  provide  a  good  description  for  ICT  excited  states  have  been  attributed  to 
the  incorrect  long-range  behavior  of  the  exchange  potential.40,41  Thus,  a  detailed  comparison  between 
computed  and  experimental  OPA  and  TPA  excitation  energies  and  intensities  for  a  series  of  coumarins 
were  carried  out  to  assess  the  accuracy  that  can  be  achieved  with  the  B3LYP,  PBEO,  and  the  new  long- 
range-corrected  (LC)  functional  of  Janai  et  al. 40,41  Previous  TDDFT  studies  on  coumarins  focused  on 
the  OPA  excitation  energy  but  not  on  intensity,  and  TPA  spectra  were  not  considered. 

The  aim  of  the  present  work  is  to  investigate  both  the  OPA  and  TPA  spectra  for  a  series  of 
coumarin  derivatives  with  different  degrees  of  ICT  characters  (Figures  1  and  2)  in  the  gas-phase  and  in 
solvents.  In  7-aminocoumarins,  the  amine  and  carbonyl  moieties  are  expected  to  act  as  electron  donor 
and  acceptor,  respectively.  Thus,  the  4-methyl  coumarin  derivatives  (Figure  1)  with  increasing  degree 
of  alkylation  measure  the  effects  of  different  donor  groups  on  the  OPA  and  TPA  spectra.  Effects  of  the 
acceptor  are  also  gauged  by  4-trifluoromethyl  coumarin  derivatives  (Figure  2). 

2.Computational  Methods 

The  computation  of  excitation  energies  and  cross  sections  was  preceded  by  the  determination  of 
equilibrium  geometries,  and  then  single  point  excited  state  calculations  were  carried  out  to  obtain  the 
vertical  excitation  energies,  transition  dipoles,  and  oscillator  strengths.  All  structures  have  been 
predicted  using  the  Kohn-Sham  (KS)42  density  functional  theory  with  the  6-31  lG(d,p)43  basis  set  and  the 
B3LYP  hybrid  functional.34"’6  These  structures  were  verified  to  be  minima  on  the  potential  energy 
surface  with  hannonic  frequency  calculations  using  the  Gaussian  03  program.44  TDDFT  calculations 
were  carried  out  using  the  6-31  lG(d,p),  6-31 1++G(d,p),  and  aug-cc-pVTZ  (for  selected  cases)  basis  sets 
with  the  B3LYP,  PBEO,  and  the  new  hybrid  functional  of  Yanai  et  al.,40  as  implemented  in  the  Dalton 
program.4 '  This  hybrid  functional,  called  CAMB3LYP,  combines  properties  of  B3LYP  and  the  long- 
range-correction  for  the  exchange  potential  by  an  Ewald  split  of  the  rn'  operator  into40 
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(1) 


_  l-fa  +  ^-erf(yur12)1  a+/Terf(/^7) 


rn  = 


rn 


where  /j,  a,  j3  are  empirical  parameters.  Yanai  et  al.  found  that  the  CAMB3LYP  functional  yields  the 
smallest  error  (MAE  =  2.53  kcal/mol)  with  (j.  =  0.33,  a  =  0.19,  and  p  =  0.46  (a  +  p  =  0.65)  in 
atomization  energies  for  the  G2  set.  This  is  slightly  better  than  the  corresponding  MAE  of  2.68  kcal/mol 
for  B3LYP,  which  is  nearly  equivalent  to  setting  a  =  0.2  and  P  =  0.0  in  CAMB3LYP  41  However,  the 
ICT  excitation  energies  for  a  model  dipeptide  are  greatly  improved,  when  compared  to  the  CASPT2 
results.  In  addition,  we  tested  a  modified  CAMB3LYP  (m C AMB 3 L YP)  functional  with  a  two- 
parameter  Ewald  split  of  the  r~2  operator  as, 

_-i  _!-[«+« -erfi!>12)]  _  a  +  a-erfQr13) 

rn  ~  +  \ 

^ i 


where  fi  and  a  are  kept  unchanged. 

To  account  for  the  solvent  effects,  TDDFT  calculations  were  carried  out  with  the 
nonequilibrium  polarizable  continuum  model46,47  (PCM)  and  the  self-consistent  reaction  field  (SCRF) 
model  of  Mikkelsen  et  al.,48  which  uses  a  spherical  solute  cavity  (SCRF-S).  The  PCM  solute  cavities 
were  carried  out  using  the  united  atom  topological  model  and  empirical  solvent  radii,  as  implemented  in 
the  Gaussian  03  program,44  The  radii  (a0)  for  the  SCRF-S  were  obtained  using  the  largest  interatomic 

distances  of  the  solutes  with  added  adjustments  for  van  der  Waals  interactions.  The  SCRF-S 
calculations  were  carried  out  using  the  Dalton  program.4" 

The  TPA  cross  sections  can  be  obtained  by  relating  the  absorption  rateJ  ’49  to  the  TPA  transition 
probability,  which  was  first  derived  by  Goppert-Mayer  using  second  order  perturbation  theory.50  The 
two  photon  absorptivity  dfo  between  the  ground-state  (0)  and  the  final  state/is  written  asJ 7,5 1,52 


Sf,{Ex,E2)-  E]  E2 g /-Q (-fi) i E2  jjs1  f0 (E1 , £2)j  » 

\ch )  / 


(3) 


where  c  is  the  speed  of  light,  h  is  the  Planck’s  constant,  fs\  the  photon  energies,  gj^EyE^  is  the 
normalized  line  shape  function  for  the  ground  to /transition,  which  can  be  represented  by  a  Gaussian  or 
a  Lorentzian.  The  two-photon  matrix  elements  SpiE^Ey)  were  obtained  from  single-  (SRQR)  or 
double-  (DRQR)  residue  of  the  quadratic  response  function  of  the  electric  dipole  operator^  within  the 
framework  of  TDDFT.54  TPA  cross  section  calculations  with  DRQR  are  carried  out  using  the  sum-over¬ 
state  (SOS)  method.  These  calculations  were  carried  out  the  the  Dalton  program. 
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larger  than  the  corresponding  MAE  obtained  with  the  PBEO  and  B3LYP  functionals.  As  illustrated  for 
Cl 52,  the  large  deviations  between  B3LYP  and  CAMB3LYP  can  be  attributed  to  the  value  of  (3  (see 
Table  IS  in  the  Supporting  Information  Section),  since  j3  =  0  for  B3LYP  (a  =  0. 19  nearly  the  same  as  in 
B3LYP).  Thus,  the  a  and  (3  parameters  must  be  reexamined  for  both  the  ground  (e.g.  the  G2  set)  and 
excited  states.  However,  the  effects  of  f3  on  the  ground  state  energetics  are  found  to  be  small.  The  MAE 
of  atomization  in  the  G2  set  for  (3  =  0  and  f3  =  0.46  changes  from  2,68  kcal/mol  to  2.53  kcal/mol, 
respectively.'10  Thus,  the  corresponding  error  for  f3  =  a  =  0.19  (Eq.  2)  is  likely  within  this  range.  The 
small  compromise  in  the  ground  state  energetics  results  in  a  significant  improvement  in  the  predicted 
Sq-Si  ICT  energies,  which  have  MAEs  of  0.10  and  0.06  eV  (without  Cl 53  due  to  convergence  problem) 
for  the  6-3 1 1  G(d,p)  and  6-3 1 1  ++G(d,p)  basis  sets,  respectively.  It  is  worth  noting  that  excitation 
energies  computed  with  the  nonhybrid  BLYP  functional  and  Tamm-Dancoff  approximation  are 
computationally  more  economical  but  do  not  offer  comparable  accuracy  delivered  by  TDDFT  with  the 
hybrid  B3LYP,  PBEO,  or  mCAMB3LYP  functional.  The  predicted30’31  gas-phase  S0-Sj  energies  of  2.99 
and  2.90  eV  for  C152A  (C35)  and  Cl 53,  respectively,  are  significantly  lower  than  the  corresponding 
experimental  values6  of  3.56  and  3.37  eV.  For  Cl 51,  the  BLYP  and  SAOP  gas  phase  excitation 
energies  of  3.32  5  and  3.24  eV,  respectively,  were  reported.  These  predicted  values  are  much  lower 
than  the  experimental  value  of  3.70  eV  in  n-hexane.12 

We  now  discuss  the  trends  observed  with  the  predicted  results  in  the  gas -phase,  which  do  not  appear 
to  be  dependent  upon  the  choice  of  basis  and  functionals.  For  example,  the  B3LYP  and  PBEO  (0.17  eV) 
S0-S j  shifts  for  C152-C153  are  identical  for  the  6-31 1G**  and  6-31 1++G**  basis  sets.  The  predicted 
value  of  0.17  eV  is  in  good  agreement  with  the  corresponding  experimental  value  of  0.22  eV.  All  four 
functionals  predict  the  Sq-Sj  energy  to  shift  to  the  red  with  increasing  degree  of  alkylation.  Cl 20  <C1 
<C480.  A  tri fluorine  substitution  at  the  4-methyl  group  gives  pronounced  red  shifts  that  increase  in 
magnitude  with  higher  degree  of  alkylation  (C151  <  C307  <C152  <152A  <  C522  <C153),  in  agreement 
with  the  experiment  data6  for  the  latter  four  coumarins.  We  should  note  that  the  absorption  maxima  of 
Cl 52  and  C152A  listed  in  Table  1  of  Ref.  6  must  be  interchanged  for  them  to  be  consistent  with  their 
corresponding  absorption  spectra  given  in  Figure  3. 

The  continuum  solvation  models  do  not  treat  specific  solvent-solute  interactions  such  H-bonding, 
which  is  important  in  water8  24  and  methanol6'  but  not  in  a  cyclohexane  solution.  For  the  6-31  lG(d,p) 


8 


basis  set,  the  SCRF-S  model  of  cyclohexane  overestimates  the  experimental  Sq-Sj  energies  by  0.05, 

0.22,  and  0.42  eV  for  the  B3LYP,  wCAMB3LYP,  and  CAMBLYP,  respectively.  Using  the  same  basis 
set,  the  errors  (MAE)  for  PCM  with  PBEO  and  B3LYP  are  0.03  and  0.1 1  eV,  respectively.  These 
excitation  energies  do  not  appear  to  be  strongly  affected  by  the  addition  of  diffuse  functions,  which  are 
found  to  lower  the  energies  by  0.05-0.07  eV  for  C120  (see  Table  1).  Thus,  the  6-311G(d,p)  basis  set 
seems  adequate  for  predicting  Sq-S  j  energies.  Note  that  the  errors  for  B3LYP  are  slightly  improved 

compared  to  those  obtained  in  the  gas-phase  due  to  the  effects  of  compensation  in  excitation  energy  by 
the  SCRF-S  method  and  the  B3LYP  functional.  This  might  be  attributed  to  the  tendency  to  overestimate 
and  underestimate  the  excitation  energy  by  the  SCRF-S  model  and  the  B3LYP  functional,  respectively. 
Consequently,  the  mCAMB3LYP  and  CAMBLYP  energies  that  are  overestimated  in  the  gas-phase,  are 
exacerbated  in  the  SCRF-S  model.  For  B3LYP,  the  SCRF-S  excitation  energies  are  larger  (by  0.03-0.07 
eV)  compared  to  the  corresponding  PCM  (with  MAE  of  0.03  eV)  values.  This  suggests  that  excitation 
energies  of  mCAMB3LYP  and  CAMBLYP  are  likely  improved  with  the  PCM  model  and  the  use  of 
diffuse  functions,  assuming  the  reduction  in  excitation  energy  observed  in  C120  holds  for  other 
coumarins.  Thus,  using  PCM  with  more  flexible  basis  sets  might  bring  the  mCAMB3LYP  errors  closer 
to  those  obtained  with  PBEO  (MAE  of  0.11  eV).  However,  it  is  unlikely  that  using  PCM  with  larger 
basis  sets  would  bring  the  CAMB3LYP  errors  to  an  acceptable  range  for  these  systems.  Overall,  we 
found  the  predicted  B3LYP  and  PBEO  Sq-Sj  energies  in  cyclohexane  are  in  excellent  agreement  with 
experiment.  The  error  for  wCAMB3LYP  and  CAMB3LYP  are  slightly  larger  than  those  obtained  in  the 
gas-phase.  The  trends  predicted  for  the  gas-phase  S0-Sj  energy  to  shift  to  the  red  with  increasing  degree 

of  alkylation  for  the  nonfluorinated  (C120  <C1  <C480)  and  for  the  trifluorinated  series  (Cl  51  <  C307 
<052  <152A  <  C522  <053)  also  hold  true  in  cyclohexane,  and  are  confirmed  by  a  collection  of 
experimental  values  in  different  studies. 

For  ethanol,  the  SCRF-S  and  the  PCM  models  appear  to  overestimate  S0-S  j  energies,  suggesting  the 
H-bonding  effects  are  significant,  especially  in  020  and  051,  The  overall  errors  for  PCM,  however. 
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remain  low  for  B3LYP  (0.1 1  eV)  and  PBEO  (0.21  eV).  The  MAEs  using  the  SCRF-S  model  are  0.13, 
0.27,  and  0.48  eV  for  the  B3LYP,  mCAMB3LYP,  and  CAMBLYP,  respectively.  The  latter  two 
functionals  overestimate  all  the  Sq-S j  excitation  energies,  with  Cl 20  having  the  largest  error.  In  fact, 
the  continuum  solvation  models  fail  to  predict  the  observed  blueshift  in  the  excitation  energy  in  water 
with  respect  to  EtOH  (see  Table  1).  This  might  be  attributed  to  the  H-bonding,  which  is  not  accounted 
for  in  the  continuum  models  of  solvation.  In  C120  (and  Cl 51),  water  and  alcoholic  solvents  can  form 
hydrogen  bonds  with  the  nitrogen  (A)  or  oxygen  (B)  at  the  amine  and  the  carbonyl  groups,  respectively, 
and/or  accepting  an  H-bond  at  the  amine  hydrogens  (C).55'57  The  B  and  C  interactions  strongly  stabilize 
the  positive  and  negative  charges  on  the  amino  nitrogen  and  carbonyl  oxygen,  respectively,  in  ICT 
excited  state  of  Cl 20,  leading  to  a  red  shift  in  the  ICT  excitation  energy  in  the  more  acidic  water 
compared  to  alcohols.  The  A  and  B  interactions  are  present  in  all  7-aminocoumarins  while  the  C 
interactions  are  unique  to  the  non-alkylated  and  the  N-alkylated  (not  N-dialkylated)  systems.  However, 
the  blueshift  in  the  S0-Sj  excitation  energy  in  water  with  respect  to  alcohols  has  been  attributed  to  the  A 
interaction  that  better  stabilizes  the  ground  state  than  the  ICT  excited  state.36,57  Interestingly,  previous 
DFT  studies  of  Cl  51  with  one  water  molecule  did  not  locate  a  stable  A-type  conformer  in  the  ground 
state.  Therefore,  calculations  of  Cl 20  in  water  and  methanol  clusters  embedded  in  a  dielectric 
continuum  are  important  to  shed  more  light  on  the  effects  on  H-bonding  on  the  solution  spectra.  Thus, 
we  have  carried  out  PCM-TDDFT  calculations  on  Cl 20  with  small  number  of  water  and  methanol 
molecules.  Although  some  conform ers  with  the  A  and  C  interactions  possess  blueshifts  in  excitation 
energies,  the  computed  Boltzmann  averaged  excitation  energies  (neglecting  entropy)  are  similar  for 
water  and  methanol.  The  inclusions  of  entropy,  which  requires  more  sophisticated  treatments  for  low 
frequencies  modes  beyond  the  so  called  rigid  rotor  and  harmonic  oscillator  approximation,  are  under 
consideration.  Calculations  for  larger  clusters  using  Monte  Carlo  optimization3*1  are  also  being  carried 
out  and  will  be  reported  in  the  future. 
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We  now  consider  the  OPA  intensity,  which  can  be  measured  with  the  Pge  or  /  values  that  are 

available  for  C152'4  and  Cl 53 11  in  a  number  of  solvents.  For  Cl 53,  the  emission  and  absorption 
transition  moments  were  reported  to  be  the  same  to  within  the  experimental  uncertainty.11  For  other 
systems,  the  pge  and  /  values  in  solutions  were  obtained  by  using  published  maxima  of  absorption 

coefficients  and  corresponding  full-width  at  half  maxima  (FWHM),  and  by  integrating  the  appropriate 
spectral  bands,  which  were  approximated  by  Gaussians.  Although  gas-phase  absorption  spectra  for  a 
number  of  coumarins  have  been  reported,  the  corresponding  absolute  intensities  were  not  measured, 
except  for  the  /value  of  0.37  for  Cl  53,'  which  is  within  the  range  (0.33-0.44)  of  the  computed  values. 
For  other  systems,  our  comparison  between  theory  and  experiment  is  largely  limited  to  values  obtained 
in  alcoholic  solutions.  Experimental  /  values  in  cyclohexane  and  water  are  also  reported  for  some 
systems. 

The  value  of / for  Cl  20  in  water  (0.325)  was  obtained  by  digitization  of  the  OPA  spectrum  reported 
by  Overway  et  al,ls  The  corresponding  value  in  ethanol  (0.318)  was  obtained  from  the  spectral  data 
reported  by  Fisher  et  al.16  and  Reynolds  and  Drexhage.60  These  experimental  oscillator  strengths  are 
slightly  smaller  than  the  computed  values  of  about  0.4-0.5,  with  those  predicted  by  CAMB3LYP  and 
mCAMB3LYP  functionals  being  the  larger  ones  (see  Table  1).  The  transition  dipoles  and  oscillator 
strengths  obtained  with  the  PCM  model  are  also  slightly  larger  than  the  corresponding  SCRF-S  values. 
The  overall  errors  for  PCM,  however,  remain  low  for  B3LYP  (0.41  D)  and  PBE0  (0.55  D).  The  MAEs 
using  the  SCRF-S  model  are  0.28,  0.48,  and  0.71  D  for  the  B3LYP,  mCAMB3LYP,  and  CAMBLYP, 
respectively.  It  is  encouraging  that  the  predicted  mCAMB3LYP  transition  dipoles  are  in  accordance 
with  those  of  B3LYP  and  PBE0.  The  predicted  changes  in  jige  from  cyclohexane  to  MeOH/EtOH  are 

small.  This  is  consistent  with  the  lack  of  dependence  on  solvents  for  absorption  transition  moments  of 
Cl 02  and  Cl  53,  as  observed  by  Lewis  and  Maroncelli.11 
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Upon  alkylation,  the  computed  |aee  and  /  values  of  the  ICT  band  show  a  slight  increase,  as  predicted 

by  the  PCM  model  using  the  PBEO  (0.07,  0.61  D)  and  B3LYP  (0.07,  0.64  D)  functionals  for  Cl.  This 
is  in  good  agreement  with  the  experimental  shift  of  0.08  and  0.74  D  for  the  oscillator  strength  and 
transition  moment,  respectively.  Note  that  these  observed  shifts  are  also  well  reproduced  using  the 
SCRF-S  model  with  the  B3LYP  (0.08,  0.68  D)  and  the  mCAMB3LYP  (0.1 0,  0.74  D)  functionals  but  not 
as  accurate  with  the  CAMB3LYP  (0.12,  0.79  D)  functional.  The  alkylation  of  Cl 20  to  create  C480 
produces  a  similar  increase  (0.09,  0.95  D)  in  the  experimental / and  pge  values  from  the  absorption  data 

of  Fisher  et  al. 16,19  and  of  Reynolds  and  Drexhage.60  However,  the  corresponding  computed  values  for 
the  syn  and  anti  conformations  of  C480  show  smaller  increases  from  that  of  Cl 20,  which  is 
overestimated  by  the  continuum  solvation  models.  Furthermore,  the  two  isoenergetic  syn  conformations 
are  0.5  and  0.2  kcal/mol  lower  in  energy  than  their  anti  counterparts  in  the  gas-phase  and  in  MeOH, 
respectively.  Thus,  rapid  syn-anti  conversions  may  occur  at  room  temperature,  leading  to  a  different 
mixture  of  conformers  and  corresponding  intensities  for  the  ICT  band.  Similarly,  the  syn  conformers  of 
the  trifluorine  substituted  Cl 53  are  0.3  kcal/mol  lower  in  energy  than  their  anti  counterparts  in  the  gas- 
phase.  This  is  identical  to  the  syn-anti  relative  energy  reported  by  Cave  et  al.29  However,  recent 
calculations  by  Sulpizi  et  al.j0jl  for  Cl 53  using  the  nonhybrid  BLYP  functional  with  a  plane-wave 
basis  set  and  pseudopotentials  revealed  that  the  anti  conformer  has  the  lowest  energy.  The  other  anti 
conformer  was  located  at  0.8  kcal/mol  above  the  lowest  anti  congener.  The  two  syn  conformers  were 
found  to  be  0.8  and  1.1  kcal/mol  above  the  lowest  anti  structure.  We  found  that  the  syn  and  anti 
conformers  of  Cl  53  are  isoenergetic  in  EtOH. 

The  shifts  in  transition  moment  for  the  trifluorine  derivatives  follow  the  alkylation  trends  discussed 
above,  increasing  from  C151  (4.97  D)  <  C307  (5.52  D),  Cl 52  (5.31  D)  <  C152A  (5.88  D),  Cl 53  (5.71 
D).  The  predicted  transition  moments  (5.23-5.69  D)  for  Cl 51  are  significantly  larger  than  experimental 
value  (4.97  D)  in  EtOH  due  to  the  importance  of  H -bonding.  For  Cl 53,  the  predicted  transition  values 
vary  from  5.84  D  (B3LYP)  to  6.40  D  (CAMB3LYP).  The  corresponding  experimental  value  of  5.71  D 
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was  obtained  by  using  the  peak  extinction  coefficient  reported  by  Reynolds  and  Drexhage611  and  the 
FWHM  reported  by  Gustavsson  et  al.12  This  value  is  in  good  agreement  with  the  transition  dipole 
moments  of  5.03-5.71  D  reported  by  Lewis  and  Maroncelli11  for  Cl 53  in  methanol,  as  obtained  by 
numerically  integrating  the  spectrum  and  scaling  them  with  different  refractive  index  correction  factors. 
Their  non-scaled  value  of  5.80  D  for  the  emissive  transition  dipole  moment,  which  was  reported  to  be 
the  same  as  the  absorptive  moment,  is  slightly  (0,1  D)  larger  than  the  value  integrated  with  Gaussian 
lineshape  above.  The  predicted  errors  of  0.3-0. 4  D  for  B3LYP  and  «rCAMB3LYP  are  close  to  the 
estimated  experimental  uncertainty  of  0.3  D  (5%)  for  Cl 53. 1 1 

B.  Two-Photon  Absorption 

The  TP  A  absorption  cross  sections  are  listed  in  Tables  2,  3,  and  4  for  various  basis  sets  and 
functionals  along  with  available  cross  sections.  In  addition,  cross  sections  obtained  with  experimental 
(in  Table  1)  and  the  difference  between  the  ground  and  first  excited  state  (Ap^e)  dipole  moments 

using  the  two-state  approximation  are  also  listed  for  four  coumarins.  These  experimental  Apge  have 

been  estimated  by  several  experimental  methods,  including  time-resolved  microwave  dielectric 
absorption  measurements,61  which  is  claimed  to  be  superior.  Note  that  these  experimental  studies  were 
not  carried  in  the  same  solvents  as  in  the  TP  A  studies.  Thus,  the  cross  sections  obtained  from  these 
quantities  offer  qualitative  rather  than  quantitative  comparisons  with  the  TPA  measurements.  The 
consensus  in  OP  A  intensity  lineshape  at  room  (or  higher)  temperature  appears  to  be  Gaussian. 
However,  both  Gaussian  and  Lorentzian  functions  are  often  used  to  compute  TPA  cross  sections.  The 
difference  between  two  line  shape  functions,  as  illustrated  for  Cl  20  in  Table  2,  is  a  factor  of  1.48,  with 
the  Gaussian  cross  (d>G)  section  being  larger  than  the  corresponding  Lorentzian  (S L)  value.  TPA  cross 
sections  computed  with  the  single  residue  of  the  quadratic  response  function  are  denoted  by  SRQR. 
DRQR-SOS  refers  to  values  computed  with  the  SOS  method  using  the  DPT  ground  state  dipoles, 
transition  dipoles  and  excitation  energies  from  linear  response  TDDFT,  excited  dipoles  and  excited-to- 

excited  transition  dipoles  from  the  double  residue  quadratic  response  function.  The  latter  approach  is 
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computationally  more  demanding,  but  offers  insights  into  the  mechanism  of  TPA.  For  Cl 20  with  the 
lowest  five  excited  states  in  the  SOS  expression,  the  two  methods  yield  cross  sections  at  the  TPA 
resonance  within  a  few  GM  of  each  other  (see  Table  2).  Our  SOS  calculations  reveal  a  type  II  TPA  for 
the  first  state  of  Cl 20,  in  which  the  TPA  is  driven  by  Apge.  Thus,  the  difference  in  the  cross  sections 
obtained  from  the  two-state  and  the  five-state  approximation  is  (1-3  GM)  small  (cf.  Tables  2  and  3).  We 
found  that  the  cross  sections  are  not  significantly  affected  by  the  basis  set.  For  C120,  changes  of  SG 
upon  going  from  the  6-31G(d)  to  the  6-31 1++G(d,p)  basis  set  using  the  SRQR  method  are  less  than  one 
GM  for  all  functionals.  For  the  largest  basis  set,  variations  in  the  peak  cross  sections  are  within  2  GM, 
as  obtained  with  CAMB3LYP  (5.5  GM),  PBEO  (6.9  GM)  ,  and  B3LYP  (7.1  GM)  functionals.  Small 
variations  among  the  functionals  are  also  observed  in  solvents  for  Cl 20  and  other  systems. 

In  solvents,  the  predicted  cross  sections  of  13-14  GM  for  Cl 20  (Table  4)  are  significantly  larger  than 
those  in  the  gas-phase,  due  to  the  increase  in  Apgc  (Table  3).  The  cross  sections  predicted  by  PCM  are 
significantly  larger  than  those  obtained  with  the  SCRF-S  using  the  two-state  approximation.  This  can  be 
attributed  to  larger  excited  state  dipole  moments  predicted  by  PCM,  It  is  not  straightforward  to  compare 
these  computed  values  with  experiment,  since  the  TPA  cross  sections  of  19  and  14  GM  measured  at 
3.29  and  4.13  eV  in  ethanol  and  water,  respectively,  are  not  the  local  maxima  of  the  TPA  spectra. 
The  maxima  are  likely  located  at  higher  (lower)  energy  with  even  larger  cross  sections  than  the  predicted 
values  for  ethanol  (water).  Interestingly,  these  cross  sections  fall  in  the  range  (2-26  GM)  of  values 
obtained  for  the  two-state  approximation  using  experimental  pge  and  Ajage. 

The  predicted  TPA  cross  section  (energy)  increases  (decreases)  with  increasing  degree  of  alkylation. 

This  is  consistent  with  the  increase  in  Ajuge  and  pge  upon  alkylation  that  has  been  predicted  and 

experimentally  observed.6'  However,  the  increase  in  the  cross-section  upon  increasing  degree  of 

alkylation  (Cl 20  <C1  <C480)  is  not  quantitatively  reproduced.  The  experimental  cross  sections  of  104 

and  160  GM  for  Cl  and  C480  (Table  4)  are  significantly  larger  than  the  corresponding  predicted  values 

of  28  (Cl)  and  22  (C480)  GM.  These  discrepancies  can  be  attributed  to  the  uncertainty  both  in  the 

14 


predicted  and  observed  cross  sections.  In  addition  to  the  errors  from  the  line  width  functions  as 
discussed  above,  the  computed  errors  might  be  attributed  to  the  specific  solvent  interactions  that  are  not 
accounted  for  in  the  SCRF  model.  Experimentally,  the  errors  for  TPA  cross  sections  for  Cl  and  Cl 20 
have  been  reported  of  about  half  order  of  magnitude.16  For  C480  (Cl 02),  the  TPA  cross  section  was 
obtained  by  relative  fluorescence  measurement19  assuming  unity  fluorescence  quantum  yield  for  both 
the  p-bis-o-m  ethyl -styrylbenzene  (BMSB)  reference  and  C480,  which  was  reported  to  be  0.5860  to  0.95 9 
in  EtOH.  Very  high  (0.995-0.870)  quantum  yields  were  reported  for  BMSB  in  a  few  nonpolar  solvents, 
but  was  not  given  for  EtOH.6'  Thus,  an  adjustment  the  TPA  cross  section  of  C480  based  on  measured 
quantum  yields  w'ould  be  useful. 

For  the  trifluorinated  series,  the  predicted  TPA  cross  sections  follow  the  alkylation  trends  previously 
observed.  The  TPA  maximum  (£G)  cross  section  (energy)  is  predicted  (SCRF-S-B3LYP/6-31 1G**)  to 
increase  (decrease)  from  14.0  GM  (3.45  eV)  to  25.1  GM  (3.17  eV)  for  Cl  51  and  C307,  respectively. 
The  TPA  spectrum  of  C151  in  EtOH  taken  from  2.95  eV  to  3.28  eV  show  a  cross  section  of  47  GM  at 
3.28  eV,16  which  is  comparable  to  the  OP  A  maximum  of  3,24  eV  in  the  same  solvent.  The 
experimental  TPA  cross  section  of  47  GM  for  Cl 51  is  more  than  twice  maximum  cross  section  of  19q  ± 
5.5  GM  obtained  in  MeOH  at  3.20  eV  for  C307,1'  even  with  the  fluorescence  quantum  yield  q  =  1. 
Thus,  the  observed  redshift  in  the  TPA  maximum  of  Cl 51  upon  alkylation  is  consistent  with  the 
theoretical  prediction.  However,  the  decrease  in  the  corresponding  cross  section  is  in  discordance  with 
theory.  It  should  be  noted  that  the  experimental  cross  section  of  47  GM  reported  for  C 1 5 1  has  a  much 
larger  uncertainty  (half  order  of  magnitude)16  than  that  of  C 3 07.  Thus,  new  experimental  studies  on  the 
TPA  cross  sections  of  Cl  51  and  other  coumarins  would  useful  in  verifying  the  predicted  results. 

4.  Summary  and  Conclusions 

We  found  that  gas-phase  excitation  energies  for  Cl 52,  C152A,  C522,  and  Cl 53  are  well  reproduced 
with  the  6-31  lG(d,p)  basis  set,  with  the  mean  absolute  errors  of  0.03,  0.08,  0.10,  and  0.32  eV,  for  PBE0, 

B3LYP,  mCAMB3LYP  (a  =  0.19  and  [3  =  0.19),  and  CAMB3LYP  (a  =  0.19  and  p  =  0.46)  respectively. 

15 


The  corresponding  set  of  0.16,  0.05,  0.06,  and  0.29  eV  are  found  for  the  6-31 1++G(d,p),  which  has 
negligible  differences  in  Sq-S  j  excitation  energies  and  oscillator  strengths  from  the  larger  aug-cc-pVTZ 
basis  set.  All  the  four  functionals  predict  S0-Sj  energy  to  shift  to  the  red  with  increasing  degree  of 
alkylation.  Trifluorination  at  the  4-methyl  group  gives  pronounced  red  shifts  that  increase  in  magnitude 
with  higher  degree  of  alkylation.  These  predictions  are  found  to  be  in  agreement  with  experiment. 

The  6-311G(d,p)  excitation  energies  predicted  by  the  SCRF-S  model  overestimate  experimental 
S0-S(  energies  by  0.05,  0.22,  and  0.42  eV  for  the  B3LYP,  mCAMB3LYP,  and  CAMBLYP, 
respectively,  in  cyclochexane.  In  MeOH/EtOH,  the  MAEs  for  excitation  energy  (transition  moment)  are 
0.13  (0.28),  0.27  (0.48),  and  0.48  (0.71)  eV  (D)  for  the  B3LYP,  mCAMB3LYP,  and  CAMBLYP, 
respectively.  The  corresponding  MAEs  for  the  PCM  model  are  0.11  (0.41)  and  0.21(0.55)  eV  (D)  for 
B3LYP  and  PBE0,  respectively.  The  larger  errors  of  excitation  energy  predicted  in  alcoholic  solvents 
are  attributed  to  the  H-bonding,  which  is  not  accounted  for  in  the  continuum  solvation  models. 

The  first  excited  state  of  coumarins  have  large  TPA  cross  sections,  which  are  driven  by  the  change  in 
dipole  moments  between  the  ground  and  the  more  polar  first  excited  state.  We  observed  large 
enhancement  in  the  TPA  cross  section  in  polar  solvents  due  to  the  increase  in  excited  state  dipole 
moments  compared  to  the  corresponding  ground  state  values.  We  found  that  the  variations  in  the  TPA 
peak  cross  sections  are  small  (1-3  GM)  among  the  various  basis  sets  and  functionals,  including  the 
mCAMB3LYP  functional  which  significantly  reduces  the  errors  in  the  TPA  and  OPA  resonances.  The 
predicted  cross  sections  are  largely  consistent  with  experiment  to  within  the  experimental  uncertainties. 
Some  large  deviations  may  be  attributed  to  the  large  uncertainties  in  the  experimental  cross  sections. 
Thus,  new  experimental  TPA  studies  for  these  coumarins  would  be  useful  in  verifying  the  predicted 
results. 


16 


Table  1.  6-31  IG(d.p)  Excitation  energy  (E,  eV),  oscillator  strength  (/),  and  transition  dipole  moment  (p,^,  D)  for  coumarin  derivatives  compared 
with  experiment.  Values  in  parentheses  and  brackets  are  computed  with  the  6-31 1++G(d,p)  and  the  aug-cc-pVTZ  basis  sets. 


Theory  Expt 


System/Method 

E(S,) 

/ 

Moi 

Energy  (poi)  \A 

Solvent 

Cl  20 

CAMB3LYP 

4.29  (4.24) 

0.373  (0.368) 

4.79  (4.79) 

PBEO 

4.09  (4.04)  [4.03] 

0.348  (0.363)  [0.355] 

4.74  (4.87)[4.82] 

B3LYP 

4.00  (3.94)  [3.93] 

0.329  (0.344)[0.338] 

4.66  (4.80)  [4.76] 

EtOH 

°SCRF-S-CAMB3LYP 

4.07  (4.00) 

0.479  (0.497) 

5.57  (5.72) 

°SCRF-S-mCAMB3LYP 

3.90  (3.83) 

0.433  (0.451) 

5,41  (5.57) 

PCM-PBEO 

3.81  (3.74) 

0.456  (0.477) 

5.61  (5.80) 

°SCRF-S-B3LYP 

3.76  (3.69) 

0.390  (0.406) 

5.23  (5.39) 

PCM-B3LYP 

3.71  (3.64) 

0.424  (0.445) 

5.49  (5.68) 

3.51^(5.54)  [0.318]c>h 

Water 

PSCRF-S-CAMB3LYP 

(3.99) 

(0.495) 

(5.72) 

PSCRF-S-mCAMB3LYP 

(3.81) 

(0.447) 

(5.56) 

PSCRF-S-B3LYP 

(3.68) 

(0.401) 

(5.36) 

PCM-B3EYP 

3.70  (3.63) 

0.418  (0.439) 

5.46  (5.65) 

3.65  (4.85)  [0.325]b 

Cyclohexane 
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Table  1.  Cont. 


4SCRF-S-CAMB3LYP 

4.19(4.13) 

0.467  (0.477) 

9SCRF-S-mCAMB3LYP 

4.02  (3.97) 

0.434  (0.446) 

3SCRF-S-B3LYP 

3.89  (3.84) 

0.401  (0.415) 

PCM-PBEO 

3.93  (3.87) 

0.472  (0.495) 

PCM-B3LYP 

3.84  (3.77) 

0.443  (0.468) 

Cl 

CAMB3LYP 

4.09 

0.495 

mCAMB3LYP 

3.89 

0.449 

PBEO 

3.82 

0.433 

B3LYP 

3.72  (3.66)  [3.65] 

0.401  (0.409)  [0.405] 

qSCRF-S-B3LYP 

3.55 

0.452 

9SCRF-S-CAMB3LYP 

3.90 

0.598 

qSCRF-S-mCAMB3LYP 

3.69 

0.529 

PCM-PBEO 

3.57 

0.529 

PCM-B3LYP 

3.47 

0.490 

rSCRF-S-CAMB3LYP  3.99  0.606 

rSCRF-S-mCAMB3LYP  3.79 


0.546 


5.42  (5.51) 

5.33  (5.44) 

5.22  (5.33) 

5.63  (5.81) 

5.52  (5.72)  3.76a 


5.65 

5.52 

5.46 

5.33  (5.43)  [5.41] 


EtOH 


5.79 

6.36 

6.15 

6.25 

6.10  3.32a  (5.63)  [0.399]^' 

Cyclohexane 

6.33 

6.17 
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Tabic  1.  Cont. 


rSCRF-S-B3LYP 

3.62 

0.488 

PCM-PBEO 

3.67 

0.548 

PCM-B3LYP 

3.57 

0.510 

C480 

CAMB3LYP  (syn) 

3.97 

0.432 

mCAMB3LYP  (.syn) 

3.76 

0.391 

B3LYP  (.syn) 

3.60 

0.347 

°SCRF-S-CAMB3LYP  {syn) 

3.78 

0.524 

°SCRF-S-mCAMB3LYP 

(syn) 

3.57 

0.462 

°SCRF-S-B3LYP  {syn) 

3.41 

0.406 

PCM-B3LYP  {anti) 

3.33 

0.447 

PCM-PBEO  (syn) 

3.44 

0.467 

PCM-B3LYP  {syn) 

3.34 

0.433 

rSCRF-S-CAMB3LYP 

3.88 

0.529 

rSCRF-S-mCAMB3LYP 

3.67 

0.475 

rSCRF-S-B3LYP 

3.51 

0.422 

PCM-PBEO 

3.55 

0.488 

5.96 


6.27 

6.13  3.5413 

5.35 

5.23 

5.04 

MeOH 

6.05 

5.84 

5.60 
5.95 
5.98 

5.85 

6.00 
5.84 
5.63 
6.02 


3.19(5.84)  [0.412]h>n 

Cyclohexane 
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Table  1.  Cont. 

PCM-B3LYP  3.45 

C151 

CAMB3LYP  4.07 

mCAMB3LYP  3.88 

B3LYP  3.73 

TDA-BLYP5  3.32 

SAOP/TZP1  3.24 

°SCRF-S-CAMB3LYP  3.78 

°SCRF-S-mCAMB3LYP  3.59 

°SCRF-S-B3LYP  3.44 

PCM-PBEO  3.48 

PCM-B3LYP  3.38 

MD-QMMM-SAOP/TZPt  3.07 

rSCRF-S-CAMB3LYP  3.94 

rSCRF-S-mCAMB3LYP  3.75 

rSCRF-S-B3LYP  3.60 


0.455 

0.382 

0.347 

.66)  [3.67]  0.312  (0.320)  [0.317] 


0.464 

0.406 

.41)  0.356(0.382) 

0.425 
0.393 


0.457 

0.410 

0.366 


5.89 


3.43d 


4.97 

4,85 

4.70  (4.80)  [4.77] 


5.69 

5.46 

5.68  (5.44) 

5.68 

5.54 


5.53 

5.37 

5.18 


EtOH 


3.24  (4.97)  [0.304]c 

3.7Qi 


hexane 


Cyclohexane 
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Table  I.  Cont. 

PCM-PBEO  3.63  0.412 

PCM-B3LYP  3.53  0.411 

C307 

CAMB3LYP  3.90  0.429 

/t?CAMB3LYP  3.70  0.389 

PBEO  3.64  0.376 

B3LYP  3.55  0.352 

°SCRF-S-CAMB3LYP  3.54  0.521 

°SCRF-S-mCAMB3LYP  3.33  0.453 

°SCRF-S-B3LYP  3.16  0.398 

PCM-PBEO  3.36  0.461 

PCM-B3LYP  3.27  0.43 1 

rSCRF-S-CAMB3LYP  3.77  0.530 

rSCRF-S-mCAMB3  LYP  3.57  0.479 

rSCRF-S-B3LYP  3.42  0.433 

PCM-PBEO  3.47  0.480 

PCM-B3LYP  3.38  0.451 


5.66 


5.54  3.56e 

5.38 

5.27 

5.22 
5.12 

MeOH 

6.23 
5.99 

5.77 
6.02 

5.89  3.15  (5.52)  [0.364]f 

Cyclohexane 

6.09 

5.94 

5.78 
6.04 

5.93  3.37d 
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Tabic  1.  Cont. 


C152 

CAMB3LYP 

3.88  (3.86) 

0.468  (0.477) 

mCAMB3LYP 

3.66  (3.62) 

0.408  (0.405) 

PBEO 

3.58  (3.51) 

0.383  (0.379) 

B3LYP 

3.48  (3.40) 

0.352  (0.348) 

qSCRF-S-CAMB3LYP 

3.46 

0.542 

4SCRF-S-mCAMB3LYP 

3.22 

0.453 

4SCRF-S-B3LYP 

3.03 

0.387 

PCM-PBEO 

3.30 

0.462 

PCM-B3LYP 

3.20 

0.427 

C152A  (C35) 

CAMB3LYP 

3.84  (3.83) 

0.496  (0.497) 

mCAMB3LYP 

3.62  (3.59) 

0.432  (0.424) 

PBEO 

3.54  (3.47) 

0.407  (0.404) 

B3LYP 

3.44  (3.35) 

0.375  (0.372) 

TDA-BLYPs 

2.99 

4SCRF-S-CAMB3LYP 

3.51 

0.593 

5.64  (5.71) 

5.42  (5.43) 

5.31  (5.34) 

5.17(5.20)  3.59S 

EtOH 

6.43 
6.09 
5.80 
6,08 

5.93  3.13b(5.31)  [0.334]m 

5.83  (5.85) 

5.61  (5.58) 

5.51  (5.54) 

5.37(5.41)  3.56S 


6.68 


EtOH 
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Tabic  1.  Cont. 


qSCRF-S-mCAMB3LYP 

4SCRF-S-B3LYP 

PCM-PBEO 

PCM-B3LYP 

rSCRF-S-CAMB3LYP 

rSCRF-S-mCAMB3LYP 

rSCRF-S-B3LYP 

PCM-PBEO 

PCM-B3LYP 

C522 

CAMB3LYP 

mCAMB3LYP 

PBEO 

B3LYP 

aSCRF-S-CAMB3LYP 

aSCRF-S-mCAMB3LYP 

aSCRF-S-B3LYP 


.28 

0.504 

.09 

0.437 

.26 

0.493 

.16 

0.457 

.70 

0.614 

.48 

0.535 

.30 

0.468 

.37 

0.510 

.27 

0.473 

.80  (3.79) 

0.466  (0.493) 

.59  (3.57) 

0.415(0.436) 

.51  (3.44) 

0.397  (0.396) 

.42  (3.34) 

0.370  (0.368) 

.49 

0.568 

.28 

0.496 

.12 

0.437 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 


6.37 


6.10 

6.32 

6.18  3.10d(5.88)  [0.408]hJ 

Cyclohexane 

6.61 

6.37 

6.11 

6.32 

6.17  3.30d 

5.68  (5.86) 

5.52  (5.67) 

5.46  (5.51) 

5.34(5.39)  3.458 

EtOH 

6.56 

6.31 

6.08 
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Table  1.  Cont. 


PCM-PBEO 

3.23 

0.497 

PCM-B3LYP 

3.14 

0.463 

bSCRF-S-CAMB3LYP 

3.66 

0.572 

bSCRF-S-mCAMB3LYP 

3.45 

0.508 

bSCRF-S-B3LYP 

3.28 

0.452 

PCM-PBEO 

3.34 

0.508 

PCM-B3LYP 

3.24 

0.475 

C153 

CAMB3LYP 

3.72 

0.444 

mCAMB3LYP 

3.49 

0.385 

TDA-BLYP5 

2.90 

PBEO 

3.41  (3.34) 

0.359  (0.360) 

B3LYP 

3.31  (3.23) 

0.331  (0.333) 

aSCRF-S-CAMB3LYP 

3.44 

0.534 

aSCRF-S-mCAMB3LYP 

3.23 

0.456 

aSCRF-S-B3LYP 

3.06 

0.396 

PCM-PBEO 

3.13 

0.454 

6.37 


6.24  3.05d 

6.42 

6.23 

6.03 

6.34 

6.21  3.25d 

5.61 

5.41 


Cyclohexane 


5.27  (5.33) 
5.14(5.21) 

6.40 

6.11 

5.84 

6.19 


3.37S  (5.5I)[0.37]k 

EtOH 
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Table  1.  Cont. 

PCM-B3LYP 

3.04 

0.422 

6.05 

2.95d(5,71)  [0.373p 

Cyclohexane 

bSCRF-S-CAMB3LYP 

3.59 

0.541 

6,30 

bSCRF-S-mCAMB3LYP 

3.37 

0.469 

6.06 

bSCRF-S-B3LYP 

3.20 

0.407 

5.79 

PCM-PBE0 

3.24 

0.468 

6.18 

PCM-B3LYP 

3.15 

0.434 

6.03 

3.15d(5.84)  [0.408] 1 

aRef.  14.  bRef.  I8.  cRef.  16 
12.  kRef. 7.  lRef.  “,mRef. 

.  dRef. 9.  eRef.  63 . 
59.  nRef.  19. 

fRogers,  J.  E. 

private  communication.  ^Experimental  gas-phase  S0-S[,  Ref.  6.  hRef.  60.  'Ref.  64.  JRef. 

°aQ  =  1 0.00  au,  s  =  24.55,  n2  =  1 .847.  Pa^  =  1 0.0  au,  e  =  78.39,  n2  =  1 .776 
9ao  =  i  1 .00  au,  8  =  24.55,  n 2  =  1 .847.  rao  =  1 1 .0  au,  8  =  2.023,  n2  =  2.028 

sRef. 31 .  lRef.  32,  excitation  energy  in  n- hexane  was  obtained  from  an  average  of  400  snapshot  of  the  MD  simulations. 
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Table  2,  S\  TPA  cross  section  {in  GM,  using  FWHM  =  0.2  eV)  and  excited  state  maximum  (in 
parentheses)  for  C 1 20a 


SRQR  DRQR-SOS  (5) 


Method 

<?G 

CAMB3LYP/6-3 1 G* 

5.1  (4.34) 

3.6(4.34) 

6.8(4.35) 

4.7  (4.35) 

PBE0/6-31G* 

6.4(4.13) 

4.4(4.13) 

6.9  (4.14) 

4.7  (4.14) 

B3LYP/6-31G* 

6.6  (4.03) 

4.6  (4.04) 

6.9  (4.05) 

4.8  (4.05) 

CAMB3LYP/6-3 1 1G*  * 

5.2  (4.30) 

3.7(4.30) 

6.5  (4.31) 

4.4(4.31) 

PBE0/6-31 1G** 

6.5  (4.10) 

4.5(4.10) 

6.5  (4.12) 

4.5  (4,11) 

B3LYP/6-31 1G** 

6.7(4.01) 

4.7  (4.01) 

6.8  (4.02) 

4.7(4.01) 

CAMB3LYP/6-31 1++G** 

5.5  (4.25) 

3.9  (4.25) 

5.9  (4.26) 

4.0  (4.26) 

PBE0/6-31 1++G** 

6.9  (4.06) 

4.8  (4.06) 

6.8  (4.07) 

4.7  (4.07) 

B3LYP/6-31 1++G** 

7.1  (3.96) 

5.0(3.96) 

7.0  (3.97) 

4.9  (3.97) 

bSCRF-S-B3LYP/6-311++G** 

13.9  (3.69) 

9.5  (3.69) 

17.4(3.70) 

11.7(3.69) 

bSCRF-S-CAMB3LYP/6-3 1 1++G** 

13.2  (4.00) 

9.0  (4.00) 

15.5  (4.01) 

10.6(4.01) 

cSCRF-S-B3LYP/6-3 1 1 ++G** 

14.1  (3.70) 

9.6(3.70) 

17.6(3.72) 

11.9(3.71) 

cSCRF-S-CAMB3LYP/6-31 1++G** 

13.8  (4.01) 

8.9(4.01) 

15.3  (4.03) 

10.5  (4.02) 

Experimental  TPA  cross  section  of  19  GM  (not  a  maximum  value)  measured  at  3.29  eV  in  ethanol,16; 
experimental  TPA  cross  section  of  14  GM  (not  a  maximum  value)  measured  at  4.13  eV  in  water.18 

bag  =  10.00  au,  s  =  78.39.  n2  —  1.776  (water) 
ca0  =  1 0  au,  e  =  24.55,  n2  =  1 .847  (EtOH) 
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Table  3.  Ground  and  excited  state  dipole  moments  (in  D)  and  TPA  cross  sections  for  Cl 20  using  the 
two-state  approximation. 


Method 

Po 

Pst 

BS2 

PS3 

Solvent. 

CAMB3LYP/6-3 1 G* 

6.3 

8.3 

7.5 

2.6 

4.8 

PBE0/6-31G* 

6.4 

9.0 

8.3 

2.4 

7.5 

B3LYP/6-31G* 

6.3 

9.1 

8.4 

2.2 

8.2 

CAMB3LYP/6-31 1G** 

6.4 

8.3 

7.4 

2.6 

5.0 

PBE0/6-3 1 1G** 

6.4 

8.9 

8.1 

2.3 

7.6 

B3LYP/6-31 1 G** 

6.4 

9.1 

8.2 

2.2 

8.3 

CAMB3LYP/6-3 1 1++G** 

6.7 

8.5 

7.4 

2.6 

4.7 

PBE0/6-3 1 1++G** 

6.6 

9.0 

8.3 

2.2 

7.3 

B3LYP/6-31 1++G** 

6.7 

9.2 

8.3 

2.1 

8.0 

Experiment3 

7.9 

EtOH 

PCM-PBEO/6-3 1 G* 

9.4 

14.5 

12.3 

3.7 

33.6 

PCM-B3LYP/6-3 1 G* 

9.4 

14.5 

12.4 

3.7 

33,0 

PCM-PBEO/6-31 1G** 

9.6 

14.7 

12.0 

11.7 

34.9 

PCM-B3LYP/6-3 1 1 G*  * 

9.6 

14.8 

12.2 

11.4 

34.4 

dSCRF-S-CAMB3LYP/6-31 1++G** 

9.5 

12.6 

10.5 

11.4 

14.5 

bSCRF-S-B3LYP/6-31 1++G** 

9.6 

13.4 

11.2 

10.5 

19.1 

Water 

PCM-B3L  YP/6-3 1 1 G*  * 

9.8 

15.2 

12.4 

11.7 

36.5 

dSCRE-S-CAMB3LYP/6-31 1++G** 

9.7 

12.8 

10.7 

11.7 

14.6 

PCM-B3LYP/6-311++G** 

10,4 

15.9 

13.0 

12.3 

42.0 

dSCRF-S-B3LYP/6-3 1 1 ++G** 

9.8 

13.6 

11.4 

10.8 

18.9 

Estimated  from  the  Lippert-Mataga  relationship  using  pg  value  of  6.0  D  and  r  value 
of  3.41  A.14 

ba0  =  10.00  au,  s  =  24.55,  n 2  =  1.847  (EtOH).  da0  =  10.00  au,  s  =  78.39,  n 2  =  1.776  (water) 
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Table  4.  Computed  TP  A  cross  section  (in  GM,  using  FWHM  =  0.2  eV)  and  excited  state  maximum  (in 
parentheses)  compared  with  experiment 


SQR 

Expt. 

System/Method 

<?G 

<5l 

5 

(TS)e 

020 

B3LYP/6-31 1++G** 

7.1  (3.96) 

5.0(3.96) 

aSCRF-S-mC  AMB3LYP/6-3 1 1 ++G** 

13.7(3.82) 

9.3(3.82) 

aSCRF -S-CAMB3L  YP/6-3 1 1 ++G*  * 

13.2  (4.00) 

9.0  (4.00) 

aSCRF-S-B3L  YP/6-3 1 1 ++G* * 

14.1  (3.70) 

9.6(3.70) 

14  (4.13)** 

4^0 

bS  CRF-S-CAMB3L  YP/6-3 11 ++G** 

13.8  (4.01) 

8.9(4.01) 

bSCRF-S-mCAMB3LYP/6-3 1 1++G** 

13.8  (3.84) 

9.4  (3.84) 

bSCRF-S-B3LYP/6-31 1++G** 

13.9  (3.69) 

9.5  (3.69) 

19  (3.29)8* 

Cl 

B3LYP/6-31 1G** 

18.6(3.73) 

12.7(3.73) 

cSCRF-S-CAMB3LYP/6-3 1 1G** 

26.2(3.91) 

17.8  (3.91) 

cSCRF-S-mCAMB3LYP/6-3 1 1G*  * 

28.0(3.70) 

19.2  (3.70) 

cSCRF-S-B3LYP/6-3 1 1G** 

28.6  (3.53) 

19.6  (3.53) 

104  (3.29)g* 

8e’k 

C480  (002) 

B3LYP/6-31 1G**  (syn) 

14.9  (3.61) 

10.2  (3.61) 

dSCRF-S-CAMB3LYP/6-31 1G**  (syn) 

21.7(3.79) 

14.9(3.79) 

dSCRF-S-mCAMB3LYP/6-31 1G**  (syn) 

22.1  (3.58) 

15.2(3.58) 

dSCRF-S-B3LYP/6-311G**  (syn) 

21.7(3.42) 

14.8(3.42) 

160  (3.17)' 

ge.l 

051 

B3LYP/6-31 1G** 

10.3  (3.74) 

7.1  (3.74) 

cSCRF-S-B3LYP/6-3 1 1G** 

14.0(3.45) 

9.5  (3.45) 

47  (3.28)S* 

14e>m 

cSCRF-S-mCAMB3LYP/6-31 1G** 

15.4(3.60) 

10.5  (3.60) 

cSCRF-S-CAMB3LYP/6-3 1 1G** 

16.8  (3.79) 

11.5  (3.79) 

C307 
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B3LYP/6-311G** 


15.4(3.56)  10.5  (3.56) 

29.0  (3.55)  19.8  (3.55) 

27.1(3.34)  18.5  (3.34) 


dSCRF-S-CAMB3LYP/6-31 1G** 


dSCRF-S-mCAMB3LYP/6-31 1G** 


dSCRF-S-B3LYP/6-31 1G** 


25.1(3.17)  17.1(3.17)  19r|  (3.20)h 


aao  =  10  au,  e  =  78.39,  n2  =  1.776  (water) 
ba0  =  1 0  au,  e  =  24.55,  n2  =  1 .847  (EtOH) 
ca0  =  1 1  au,  s  =  24.55,  n2  =  1 .847  (EtOH) 
da0  =  1 1  au,  e  =  32.63,  n 2  =  1 .758  (MeOH) 

eTPA  cross  sections  are  obtained  with  experimental  \ige  in  Table  1,  A\ige  from  time-resolved  microwave 
dielectric  absorption  measurements  in  benzene,61  computed  angles  (Cl 20  =  160°,  Cl  =  165  °,  C480  = 
164°,  C151  -  164°)  between  \x.ge  and  Ap™,  and  Gaussian  lineshapes  (FWHM  =  0.2  eV)  using  the  two- 
state  approximation. 

^Experimental  TPA  cross  section  measured  in  water. 1S 

^Experimental  TPA  cross  section  measured  in  EtOH  (relative  error  of  about  1/2  of  magnitude).16 

b  Experimental  TPA  maximum  cross  section  measured  in  MeOH,  q  is  the  quantum  yield.1'  q  =  0.61,57 
0.5660  in  EtOH.  q  =  0.90  in  50%  EtOH.9 

■Experimental  TPA  maximum  cross  section  measured  in  MeOH.19 

JA \ige  =  2.4  D,  using  Ap^e  =  4.1 -6.2  D  (electrooptic)61  gives  Sa  (TS)  =  11-26  GM  using  A yige  =  1 .90  D14 
(solvatochromism)  gives  SG  (TS)  =  2  GM.  using  Ajuge  =  3.97  D65  (solvatochromism)  gives  SG  (TS)  -  1 1 
GM. 

kAgge  =  2.9  D,  using  Afige  =  7,3-7. 8  D61  (electrooptic)  gives  SG  (TS)  =  50-57  GM,  using  Apge  =  2.72 
D15  (solvatochromism)  gives  5Q  (TS)  =  7  GM. 

*A|ige  =  3.0  D,  using  A p?(?  =  2.11  D66  (solvatochromism)  gives  SQ  (TS)  =  4  GM 

mApge  =  4.4  D,  using  A \ige  =1.7  D66  (solvatochromism)  gives  SG  (TS)  =  2  GM,  using  Apge  =  2.81  D63 
(solvatochromism)  gives  SG  (TS)  =  6  GM,  using  Apge  =  4.35-7.01  D5  (electrooptic)  gives  SQ  (TS)  =  14- 
35  GM, 

*Not  a  maximum  value 
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Cl  20 


Cl  (C481) 


C480  (syn) 


C480  (anti) 


Figure  1.  Structure  of  coumarin  120,  1,  480  (102). 
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C152 


C152A  (C35) 


C522 


C153 


C307 


C151 


Figure  2.  Structure  of  coumarin  152,  152A,  522,  153,  307,  and  151. 
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